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Abstract. It is shown that a vortex can be continuously created in a Bose-Einstein 
condensate with hyperfine spin F = 2 in a Ioffe-Pritchard trap by reversing the axial 
magnetic field adiabatically. It may be speculated that the condensate cannot be 
confined in the trap since the weak-field seeking state makes transitions to the neutral 
and the strong-field seeking states due to the degeneracy of these states along the 
vortex axis when the axial field vanishes. We have solved the Gross-Pitaevskii equation 
numerically with given external magnetic fields to show that this is not the case. It is 
shown that a considerable fraction of the condensate remains in the trap even when 
the axial field is reversed rather slowly. This scenario is also analysed in the presence 
of an optical plug along the vortex axis. Then the condensate remains within the 
F z — 2 manifold, with respect to the local magnetic field, throughout the formation of 
a vortex and hence the loss of atoms does not take place. 
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1. Introduction 

It has been observed that the Bose-Einstein condensate (BEC) of alkali atom gas 
becomes superfluid Superfluidity of this system is different from the previously 

known superfluid 4 He in many aspects. For example, the BEC is a weak-coupling gas 
for which the Gross-Pitaevskii equation is applicable while superfluid 4 He is a strong- 
coupling system. One of the most remarkable differences is that the BEC has spin 
degree of freedom originating from the hyperfine spin of the atom, and that this degree 
of freedom couples to external magnetic fields. Accordingly the order parameter of the 
condensate is also controlled at will by external magnetic fields. Superfluid 3 He also 
has similar internal degrees of freedom, which, however, are rather difficult to control 
by external fields [|J. 

Taking advantage of this observation, we proposed a simple method to create a 
vortex in a BEC with the hyperfine spin F — 1 j|, [5] ; a vortex-free BEC is intertwined 
topologically by manipulating the magnetic fields in the Ioffe-Pritchard trap to form a 
vortex with the winding number 2. This is achieved by reversing the axial magnetic 
field adiabatically while the planar quadrupole field is kept fixed. 

In the present paper, a similar scenario is analysed for a BEC with F = 2, taking 
87 Rb as an example. The difference between the present case and that for F = 1 
will be emphasised in our analysis. In the next section, we briefly review the order 
parameter of BEC with hyperfine spin F = 2 and the Gross-Pitaevskii (GP) equation 
which describes the time-evolution of the order parameter. In section 3, the ground 
state order parameter of the BEC in the weak-field seeking state confined in a harmonic 
potential is obtained. Then the time-evolution of the condensate, as the axial field is 
adiabatically reversed, is studied by solving the GP equation numerically. Cases with 
different reverse time are analysed to find the best possible reversing time for which the 
the fraction of the remaining condensate in the vortex state is maximised. It is shown 
that the condensate in the end of this scenario has the winding number 4. In section 4, 
the GP equation is solved in the presence of an optical plug along the vortex axis. The 
BEC remains in the F z = 2 weak-field seeking state, with respect to the local magnetic 
field, throughout the development, and hence no atoms will be lost during the formation 
of a vortex. Section 5 is devoted to conclusions and discussions. 

2. Order Parameter of F = 2 BEC 

2.1. General F = 2 condensate and Gross-Pitaevskii equation 

Suppose a uniform magnetic field B parallel to the z-axis is applied to a BEC of alkali 
atoms with the hyperfine spin F = 2. Then the hyperfine spin state of the atom is 
quantised along this axis; the eigenvalue m of F z takes a value —2 < m < +2, where 
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F z \m) = m\m). Let us introduce the following conventions 
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The order parameter is expanded in terms of \m) as 

2 

|*)= ^rn\m) = (* 2 , ^1,^0,^-1, ^- 2 ) T , 

m=— 2 

where T denotes the transpose. 

The representation of the angular momentum operators F k (k 
is easily obtained from the well-known formulae 

(2,m\F + \2,m') = a/(2 - m)(3 + m)5 m , m '+i, 

(2,ra|F_|2,m') = a/(2 + m)(3 - m)5 m , m /_i, 

(2,m|F z |2,m') = rn<!> m , m >, 

where F ± = ± 

The dynamics of the condensate in the limit of zero temperature is given, within 
the mean field approximation, by the time-dependent Gross-Pit aevskii (GP) equation 
with spin degrees of freedom. This equation, obtained by Ciobanu, Yip and Ho (see 
also 0), is written in components \l/ m as 
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where summations over k = x,y, z and —2 < n,p,q < 2 are understood. Here, M is the 
mass of the atom and V(r) is the possible external potential. The Larmor frequency is 
defined as ftwik = 7/x-Bfe, where \lb is the gyromagnetic ratio of the atom and \xb 
the Bohr magneton. The interaction parameters are expressed in terms of the s-wave 
scattering length ap, F being the total hyperfine spin of the two-body scattering state, 
and are given by 
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where a = 4.73nm, a<i = 5.00nm and 04 = 5.61nm for 87 Rb atoms. This should be 
compared with F = 1 BEC where there are only two types of scattering state and hence 
two interaction terms in the GP equation. 
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2.2. Weak- field seeking state 



Suppose a strong magnetic field B is applied along the 2-axis. Then the components 
with F z = 1 and 2 are in the weak-field seeking state (WFSS) while those with F z = — 1 
and —2 are in the strong- field seeking state (SFSS). The presence of two WFSSs leads 
to an interesting two-component vortex that is not observed in F = 1 BEC as we see in 
the next section. The energy of the state with F z = is independent of the magnetic 
field and will be called the neutral state (NS) hereafter. Suppose a uniform condensate 
is in the state with F z = 2. The order parameter of the condensate takes the form 

|* ) = /o(l,0,0,0,0) T (4) 

where |/o| 2 is the number density of the condensate. Now let us consider a state which 
is quantised along an arbitrary local magnetic field 
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Let Fb = B-F/Bbe the projection of the hyperfine spin vector along the local magnetic 
field. The WFSS |^) which satisfies F B \^) = +2|^) is obtained by rotating |^ ) by 
Euler angles a, j3 and 7 and is given by 
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The GP equation restricted within the WFSS is obtained by substituting Eq. 
into Eq. (§), see the next section. 



3. Vortex formation without optical plug 

The formation of a vortex in the F = 2 condensate is analysed in this and the next 
sections. In the present section, we study the scenario without an optical plug along 
the vortex axis. Although some fraction of the condensate is lost from the trap in this 
scenario, the experimental setup will be much easier without introducing an optical plug. 
In fact, it will be shown below that a considerable amount of the condensate remains in 
the trap by properly choosing the time-dependence of the magnetic field. 



3.1. Magnetic fields 

Suppose a condensate is confined in a Ioffe-Pr it chard trap. It is assumed that the trap 
is translationally invariant along the z direction and rotationally invariant around the 
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z-axis. The quadrupole magnetic field of the trap takes the form 
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where <fi is the polar angle. The magnitude B±(r) is proportional to the radial distance 
r near the origin; B±(r) ~ B' ± r. The uniform time-dependent field 
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is also applied along the z-axis to prevent Majorana flips from taking place at r 
where ,B_|_ vanishes. Now the total magnetic field is given by 



B(r,t) = B ± (r) + B z (t) 



Comparing this equation with Eq. 
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It has been shown in the previous works for F 
condensate in the beginning will end up with a condensate with a vortex of the winding 
number 2 if B z reverses its direction while B± is kept unchanged @, [5] . We expect that 
the same magnetic field manipulation to lead to the vortex formation in a BEC with 
F = 2. The uniform axial field B z (t) must reverse its direction as 

B.(0)(l-f) 0<t<T 

-B g (0) T<t. 1 ' 

to create a vortex along the z-axis. This gives a "twist" to the condensate leading to 
the formation of a vortex with winding number 4, see below. 

Before we start the detailed analysis, it will be useful to outline the idea underlying 
our scenario. Suppose one has WFSS with Fb = 2 in the trap at t = 0. The magnetic 
field at r ~ points +z direction (i.e., (3 ~ 0) and hence the WFSS takes the form \l/o 
of Eq. (|4]). Then the angle 7 must satisfy 7 = —a for the BEC to be vortex- free, see 
Eq. (||). The field B z (t) vanishes at t = T/2, for which (3 = tt/2, and the hyperfine spin 
is parallel to the quadrupole field B±. Accordingly one must choose a = —(f> for this 
condition to be satisfied, see Eqs. (0) and (p~0|) . This also implies 7 = +0. When the 
field B z is completely reversed at t = T, the magnetic field at r ~ points down and 
hence (3 ~ n there. Substituting a = —7 = <p and (3 = n into Eq. (§), one finds the 
order parameter at t — T; 

m = f e- 2l(t ' (0,0,0,0,e- 2l *) T , 

which shows that a vortex with the winding number 4 has been created. 



(12) 
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3.2. Initial state 



Suppose a vortex-free BEC is confined in a Ioff-Pritchard trap, whose magnetic field 
takes the form @ and that the condensate is in the eigenstate Fb = 2 with respect 
to the local magnetic field B with B z = B z (t = 0). The condensate wave function is 
then obtained by solving the stationary state Gross-Pitaevskii equation. Substitution 
of Eq. @ with a = —7 = — into Eq. (||D yields 

h 2 

-^V 2 (/ ^m) + (91 + 4g 2 )foV m + hu L f v m = fif V m , 

where we have put \I/ m = / u m . Note that the g 3 term vanishes identically for the present 
state. The condensate wave amplitude fo(r) is taken to be a real function without loss of 
generality. The eigenvalue \x is identified with the chemical potential. If one multiplies 
the above equation by {v m Y from the left and uses the identity J2 m \ v m\ 2 — 1 an d other 
identities derived from this, one obtains the reduced GP equation for /o( r ); 
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comes from the rotation of the five- dimensional local orthonormal frame that defines 
the order parameter. The reduced GP equation looks similar to the ordinary scalar GP 
equation except that there is an extra term f3' 2 in v^S7 2 v m . 

It is convenient to introduce the energy scale fhu and the length scale ano define by 
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For a typical values B z (0) 
1.69 • 10~ 24 erg and ano — 0.68/im. After scaling all the physical quantities by these 
units, one obtains the dimensionless form of the reduced GP equation; 



IG,B' ± = 200G/cm for 87 Rb, one obtains fku 



^'+4 + 

r 



4(3cos 2 /3-5) sin 2 ^-/3 
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(16) 



where the dimensionless quantities are denoted with tilde. For example, f = r/ano, fo — 
fo a uo an d gk = 9k/ (^ a Ho)- The tilde will be dropped hereafter unless otherwise stated 
explicitly. 

Figure 1 shows the ground state condensate wave function obtained by solving 
Eq. (16) numerically. We have chosen fo(r = 0) = 6 which yields the central density 
n ~ 1.17- 10 14 cm~ 3 . This is roughly of the same order as that realised experimentally. 
The difference between the chemical potential and the Larmor energy at the origin is 
S/j, = fi — lol = 3.95, which amounts to 5fi = 6.68 • 10 _24 erg in dimensionful form. 
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Figure 1. The initial condensate wave function fa(r) in dimensionless form. The 
radial coordinate r is also dimensionless. 

3.3. Time development 

Now the time-dependent GP equation (Q) is solved numerically with the initial condition 
= fo(r)v m with f Q been obtained in the previous subsection. We have introduced a 
tanh-shaped cutoff to mimic the loss of atoms from the trap; particles reach at L 3> ano 
vanish from the system. We have made several choices of the reversing time T and 
maximised the fraction of the condensate left in the trap in the final equilibrium state. 
The details of the algorithm are given in || and will not be repeated here. 

Figure 2 shows the wave functions \^f m \ f°r the choice T/r = 1000, where r = 2-k/ujl 
is the time scale set by the Larmor frequency at t — and r = 0. One obtains 
r ~ 7.14 • 10~ 7 sec for the parameters given in the previous subsection. The parameter 
r is expected to be the measure of the adiabaticity. There are two weak-field seeking 
states possible for F = 2, those with Fb = 2 and Fb — 1. It turns out that the 
final vortex state is a mixture of these two states. When the axial field B z (t) vanishes 
at if: = T/2, the gaps among WFSSs, SFSSs and NS disappear at r = and the 
level crossing takes place there. Then the adiabatic assumption breaks down and 
some fraction of the condensate transforms into SFSSs and NS as well as Fb = 1 
WFSS. Those components in SFSSs and NS eventually leave the trap and the final 
condensate is made of Fb = 2 and Fb = 1 components. It is a remarkable feature of 
the F = 2 BEC, compared to its F = 1 counterpart, that the vortex state thus created 
is mixture of these two WFSSs. The composite nature of the final vortex state is best 
revealed by projecting \^(r)) to Fb = 1 and 2 states. Let \v) be the vector defined 
in Eq. (§) and \u) = exp(-iaF z ) exp(-(3F y ) exp(-ijF z )\l) . Then U 2 (r) = (v(r)|*(r)) 
and IT^r) = (u(r)\^(r)) depict the projected amplitudes of \^?(r)) to the local F B = 2 
and Fb = 1 state, respectively. These amplitudes are shown in Fig. 3 for \^(r)) at 
t = 10T. It is interesting to note that the Fg = 2 component has a winding number 4 
while Fb = 1 has 3. 

The fraction of the condensate left in the trap at time t has been plotted in Fig. 4 
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Figure 2. Time dependence of the order parameter |* m | for the reversing time 
T/t = 1000. 




Figure 3. The projected amplitudes |n 2 (r)| and |IIi(r)| obtained from the order 
parameter |*(r)) at t = 10T. 
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Figure 5. The fraction of the condensate left in the trap, as a function of T/t, when 
the BEC reaches equilibrium at t 3> T. The curve is shown for guide. 

for T/t = 1000. It should be noted that ~ 2/5 of the condensate is left in the trap 
when the system reaches an equilibrium at t ^> T. 

Figure 5 shows the fraction of the condensate left in the trap in the equilibrium 
state at t 3> T for various T. It can be seen from this figure that a considerable amount 
of the condensate is left in the trap for a wide variety of the reversing time T. 

In the next section, we analyse the creation of a vortex in the presence of an optical 
plug along the centre of the system. It will be shown that the vortex thus created is 
purely made oi F B = 2 WFSS. 

4. Vortex formation with optical plug 

The loss of the condensate in the previous section takes place since the energy gaps 
among WFSSs, NS and SFSSs disappear at r = when B z vanishes at t — T/2. One 
may introduce an optical plug along the vortex axis to prevent the condensate from 
entering this "dangerous" region. An optical plug may be simulated by a repulsive 
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Figure 6. Time dependence of the condensate amplitude /o in the presence of the 
optical plug. The reversing time is T/t = 10000. The condensate amplitudes at t = 
and t = T are almost degenerate. 



potential 



V(r) = V exp - 



(17) 



where Vo is determined by the power of the blue-detuned laser while ro by its waste size. 
We take V = 9.27 • 10~ 21 erg and r = 5/im in our computation below. 
Now the time-independent GP equation is given by 

2 .. . . . ,3 
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;2 -(3 cos 2 (3-5) sin 2 £ - (3' 2 + V(r) 



fo 



(18) 



+ (^i + 4^)/ 3 + ^ L / = ^/o 

in dimensionless form, where V(r) = V{r)/fko. The angle (3 is given by (3(r) = 
tan _1 [i?^(r)/i? 2 (0)]. We will drop tilde from dimensionless quantities hereafter unless 
otherwise stated. The ground state condensate wave function is obtained by solving 
this equation numerically. We find the relative chemical potential 5/i — /i — ujl — 173, 
which amounts to 8/i = 2.92 • 10~ 22 erg in dimensionful form, and the condensate wave 
function f shown in Fig. 5. 

The time-dependent GP equation 



.dfo 



dt 



r 



^(3 cos 2 (3-5) sin 2 ^ - (3 12 + V(r) 



fo 



+ (g 1+ Ag 2 )f +uj L fo (19) 

is solved with the initial wave function obtained above. Here (3 = (3(r,t) = 
ta,n~ 1 [B±(r)/B z (r)]. It should be stressed again that the condensate remains within 
the Fb = +2 WFSS throughout the temporal evolution. Figure 7 shows the time- 
dependence of the components ty m for the choice T/t = 10000, namely T = 7.14ms in 
dimensionful form. In contrast with the case without optical plug, the time-dependence 
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of the order parameter is independent of the choice of T so far as T/r is large enough 
so that the adiabaticity is observed. 

The vortex thus obtained has a region near the origin (r ~ 0) where the condensate 
cannot approach due to the presence of the optical plug. The vortex current flows around 
a multiply connected region. This situation is analogous to the superconducting current 
flowing around a ring. It is natural to expect that a vortex without the optical plug may 
be obtained if one withdraws the optical plug after the persistent current is established 
at t — T. (Note that the optical plug has been introduced to prevent Majorana flips at 
r ~ at t ~ T/2. Accordingly the optical plug is not required anymore for t > T.) Let 
us suppose that the optical plug is slowly turned off after t = T with the time constant 

V(rt) = i VoeM~rVrl) 0<t<T 

Klr ' tJ \ V exp(-r 2 /r 2 ) exp[-(t - T)/t ] T<t. ^ 

It is found that the condensate oscillates back and forth for small to- F° r sufficiently 
large t , however, the condensate smoothly rearranges itself to a vortex state without 
the optical plug. Figure 7 shows our numerical result for T/r = to/r = 10000, for which 
one still observes such oscillations. 

A vortex with a winding number 4 is thus created without losing any atoms from the 
trap. It should be noted, however, that it is technically difficult, albeit not impossible 
||, to introduce an optical plug with a few microns of radius along the centre of the 
BEC whose radial dimension without optical plug is of the order of a few microns. 

5. Conclusions and Discussions 



The formation of a vortex in a BEC with F = 2 in a Ioffe-Pritchard trap has been 
considered by fully utilising the spinor degrees of freedom. It was shown that a vortex 
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5 10 r 15 20 

Figure 8. Time dependence of the condensate amplitude /o for the reversing time 
T/t = 10000. The potential decays exponentially with the time constant to = T for 
t > T , see Eq. (20). 

with winding number 4 is created continuously from a condensate without a vortex, by 
simply reversing the axial magnetic field B z (t). This scenario has been studied with and 
without an optical plug at the centre of the vortex. Some amount of the BEC is lost 
from the trap in the absence of an optical plug while no atoms are lost in the presence 
of it. Our numerical analysis shows that there remains a considerable fraction of BEC 
even without the optical plug. The introduction of an optical plug in a trapped BEC is 
difficult, albeit not impossible. 

Our vortex has a large winding number 4 and is expected to be unstable against 
decay into four singly quantised vortices in the absence of an optical plug. Whether 
a vortex with such a large winding number may be observable depends on how large 
the lifetime of the metastable state is compared to the trapping time of the BEC. Our 
prelimiary analysis of the Bogoliubov equation suggests that the lifetime is of the order 
of 100ms and these highly-quantised vortices exist for a considerable duration of time. 

We would like to thank Kazushige Machida and Takeshi Mizushima for discussions. 
One of the authors (MN) thanks Takuya Hirano, Ed Hinds and Malcolm Boshier for 
discussions. He also thanks Martti M. Salomaa for support and warm hospitality in the 
Materials Physics Laboratory at Helsinki University of Technology, Finland. MN's work 
is partially supported by Grand-in-Aid from Ministry of Education, Culture, Sports, 
Science and Technology, Japan (Project Nos. 11640361 and 13135215). 

Note Added - - After we submitted our manuscript, the MIT group reported the 
formation of vortices according to the present scenario f|. They employed hyperfine 
spin states F — 1 and F = 2 of 85 Rb and found that the vortex thus created had the 
winding number two in the former case while four in the latter case, in consistent with 
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our prediction. The vortex state has considerably long lifetime, at least 30ms after its 
formation, in spite of higher winding number, which suggests that these vortices are 
rather stable. The stability analysis of highly-quantised vortices is outside the scope of 
the present work and will be published elsewhere. 

We are grateful to Aaron Leanhardt for informing us of their result. 
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